Coupling of intrinsic Josephson oscillations in layered 
superconductors by charge fluctuations 
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ABSTRACT 



The coupling of Josephson oscillations in layered superconductors is studied with help of a tunneling Hamiltonian 
formalism. The general form of the current density across the barriers between the superconducting layers is derived. 
The induced charge fluctuations on the superconducting layers lead to a coupling of the Josephson oscillations in 
different junctions. A simplified set of equations is then used to study the non-linear dynamics of the system. In 



particular the influence of the coupling on the current-voltage characteristics is investigated and upper limits for the 
coupling strength are estimated from a comparison with experiments on cuprate superconductors. 
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1. INTRODUCTION 

The superconducting properties of the highly anisotropic cuprate-superconductors Tl2Ba2Ca2Cu30io+5 (TBCCO) 
and Bi2Sr2CaCu208+,5 (BSCCO) are well described by a stack of Josephson junctions coupling the superconducting 



CuC>2 layers in c-d^ciikm-In particular the multiple branch structure observed in the current-voltage characteristics 
by several groupsEJBBffiHLl can be explained by this model. 

^ . Due to the low value of the critical current in c-direction the system has a small Josephson plasma frequency to p 

' and a large value of the McCumber parameter j3 c , which causes a strong hysteretic behaviour of the current-voltage 
characteristics. The low value of the Josephson plasma frequency manifests itself in the transparency of the stack 
with respect to THz^radiation in c-direction. a The longitudinal and transversal plasma oscillations have also been 
, observed directly BE3 

The theoretical investigation of longitudinal plasma oscillations has become popular since Koyama and Tachiki0 
proposed a coupling of Josephson oscillations in different barriers due to charge fluctuations. In systems with weakly 
coupled superconducting layers the charges on different layers need not to be constant, which is iiijcpntrast to ordinary 
superconductors where charge neutrality can be assumed. In the theory of Koyama and Tachikilid it follows that the 
I . gauge-invariant scalar potential 
13 ' W = *J - (ft/2e)*, (1) 

does not vanish. Here <&i is the electric scalar potential and xi is t ne phase of the superconducting order parameter 
O ■ A; = I A; I exp(ixi) on layer I. On the other hand the Josephson current density j c sin 7^2+1 between layers I and I + 1 
depends on the gauge-invariant phase difference 

• i-H , 

^ ' 7J,l+i(*) =»(*)- Xl+iW - T / dzA z (z,t) (2) 



h ji 

where A is the vector potential (in our notation e = |e|). Its time derivative (the second Josephson relation) 

h f l+1 

2^7M+i(*) = J i dz E z + fii +1 - fii (3) 

then not only depends on the voltage between the layers, but also on the potential difference — /if. This finally 
leads to a coupling between Josephson oscillations in different barriers. A non- vanishing generalized scalar potential 
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which is related to quasi-particle charge imbalance can also be obtained ssrithout breaking charge neutrality. The 
importance of this effect has been stressed by Artcmenko and KobelkovEj and by RyndykE31i-!l In the layered 
cuprate-superconductors probably both effects are present. 

In this communication we want to derive the coupling effect in a microscopic model starting from a tunneling 
Hamiltonian. We arrive at an expression for the current between different layers which is formally similar to that 
obtained by Artemenko for a different model. In leading order in the interlayer hopping t± our results are also similar 
to the model of Koyama and Tachiki. This simplified model will be used to study the non-linear dynamics of the 
system and to discuss implications of the coupling on the current-voltage characteristics. 



2. GENERAL OUTLINE OF THE THEORY 

We start from a model where the current between superconducting layers across the insulating barrier is described 
by a time-dependent tunneling Hamiltonian 

H T = £ T k , k/ 4 +hk ,^ k ^r ldzA ^+h.c. (4) 

which depends on the vector potential A z (t) in the barrier. Here TJ.^/ is a tunneling matrix element describing 
(random) hopping between neighboring layers. In order to get a Josephson current also for d-wave superconductors 
we have to keep some angular dependence in this matrix element. 

In addition to this the current is driven by the difference of scalar potentials &i(t) on the different layers. Thus 
the total time-dependent Hamiltonian is 

H = J2(Hi-eMt)N l ) + H T (t) (5) 
l 

where Hi is the Hamiltonian of the electrons in layer I including superconducting interactions and pairing. 

Finally we have to take into account phase fluctuations of the superconducting order parameter induced by 
charge fluctuations. It can be shown that the results of the calculation for physical quantities depend only on the 
gauge-invariant combinations 7;, 2+1 and /i/ of the electromagnetic potentials with the phase of the order parameter. 
Formally the same results are obtained if we replace in the Hamiltonian the exponential by exp(i7y+i(£)), the scalar 
potential $z(t) by ni(t) and assume the order parameter to be real in the BCS-treatment of Hi- 

The non-linear Josephson effect results from the periodic dependence on the phases 7(i). In the limit of large 
McCumber parameter j3 c the phase can be written as 

7 (0 = 70 + ujt + 5j(t) (6) 

where 70 is the constant phase determined by the dc-current and u> is the Josephson frequency which is related to 
the dc- voltage. The part Sj(t), which oscillates with the same frequency u>, is small for large p c . 

We calculate the current response Jn+i between neighboring superconducting layers with respect to both the 
tunneling Hamiltonian Hj- and the generalized scalar potential /i; restricting ourselves to second order tunneling 
processes and linear response with respect to /i;. 

The result can be written in the following general form: 

JM+i(*)=i^+i(*)+iM+iW (7) 



where 

iuVi(*)=f dh S (t - cos 7(t) " lih) + f dhf dt 2 S 1 (t, h , t 2 ) cos 7 W ~ 7(tl) (t 2 ) - m (fa)) 

(8) 
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Figure 1. Graphs for the current density. Symbols: left X = current operator, right X = Ht, • = density 
vertex. Each cross corresponding to a hopping T^k' between layers I and I + 1 is combined with a phase factor 
exp(±i7 M+ i(i)/2). 
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Figure 2. Graphs for the density response in layer I. Symbols: X = Ht, • = density vertex. 



is the quasi-particle current density and 



= icsin 7 (t)+/ dh J (i-i 1 )lM C os 7W+7(<l) 



2 2 
* /"* 7(*)+7(*i) 

rfti / dt 2 Ji(t,t 1 ,t 2 ) cos {m + i(t 2 ) - m{t 2 )) 



(9) 



is the Josephson current carried by the condensate ("f(t) := 7;.;+i(i)). The functions So(t— 1\) and Jo(t— ti) result from 
a folding of two normal and anomalous Green's functions in neighboring layers (see Fig. |l|), and j c = Jo(t, t) — const 
is the critical current density. In the terms Si(t, h,t 2 ) and J\{t, t\,t 2 ) the additional linear dependence of the Green's 
functions on fJ-i(t 2 ) is considered. A similar (but more complicated expression) is obtained for the charge jdfinsity 
response (see Fig. |J). Equations with a similar structure have been obtained by Artemenko and KobelkovLLj for a 
different model. 

For a systematic evaluation one has to insert the ansatz for 7; i ;+i(i) into these expressions and has to separate 
different harmonics in the Josephson frequency lo. A considerable simplification is obtained if one keeps the non- 
linear effects only in the sin-term of the Josephson current and linearizes the other terms with respect to 7. Then 
one obtains for the current density across the barrier with thickness b 



jui+i = j c sm7 M+ i + cr o^7(,i+i + o-i (Mi - Vi+i)/b 



(10) 
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and for the density response 



s Pi = X P "Jpi + xf)(m+\ + W-i - 2w) + 2-7<ri(7i_i,i - 7y+i) • ( n ) 

In the normal state o~o = o~\ and jij+i depends only on the voltage across the barrier. In the superconducting state 
7z,z+i and pi have separate physical meaning and we need the second equation for Spi to determine pi as a function of 
the voltage. The density reponse function Xpp ls finite in the superconducting state. At low temperatures it is only 
weakly frequency dependent. Approximately it is given by Xpp — — 2e 2 A 2 (0), where N 2 (0) is the two-dimensional 

(2) 

density of states of the electron gas in the CuC>2-layers. The conductivities o"o,i as well as Xpp describe the charge 
exchange with the neighboring layers and are proportional to t\ . Adding finally the displacement current we obtain 
a relation between the electronic current density across the barrier with the external current density, 

3 = 3l,i+i + eeo-E/.i+i • ( 12 ) 



In a first step we eliminate the scalar potential difference pi in favour of the gauge-invariant phase difference 
7u + i and the electric field by using the Josephson relation 

^7i,i+i(<) = bE u+ i + pi+i - pi . (13) 
In the next step we express the charge fluctuations with help of the Maxwell equation 

Spi = e e(Eij +1 - Et-xj) . (14) 
We then finally arrive at the following differential equation for the phase: 

j = (l - «A< 2 ') sin Ti.i+1 + i (l - »/A< 2 >) 7u+i + ^ (l - CA (2) ) 7 W +i ( 15 ) 

where w 2 — j c 2eb / '{ee^h) , 1/uj c = ero/(eeoWp) The dimensionless quantities a, 77, Q describe the coupling of the phase- 
difference in different layers via the derivative operator A^ 2 ), which is defined as fi = fi + \ + — 2fi. In 
particular 

a = -ee /(bx° pp ) + O(tl), (16) 
r, = -eeo/(b X pp )(l - 2a x /a Q ) + 0(t 2 ± ) . (17) 

C is proportional to t\. For to -C A the quantity a is only weakly frequency dependent (see Fig. |^). Therefore it 
will be approximated by its value at id = 0, a(0) = eeo/(2e 2 WV 2 (0)). If we neglect rj and £ we arrive back at the 
theory of Koyama and Tachiki.til In fact at uj <C T <C A we find 01 ~ ctq/2 for d-wave superconductors. Thus 
neglecting rj and ( seems to be a good approximation for small values of u>. For strong coupling in c-direction one 
ends up with a situation, where charge ihictuations on the layers are suppressed, Spi — 0, but a finite pi is generated 
by charge-imbalance of quasi-particles.Eil A study of these effects in this model will be done in the future. 



3. INFLUENCE ON THE CURRENT- VOLTAGE CHARACTERISTICS 

In our further discussion of the experimental implications of this effect on the current-voltage characteristics we 
restrict ourselves to the approximation (|T^) with a — const, 77 = £ = supplemented by 

^-7i,i+i = Vu+i - a(Vi +lii+2 + Vj_i,j - 2V hi+1 ) (18) 
ze 

where we have defined the voltage Vy+i = f, + E z dz = bE^i+i . In particular, we obtain a relation for the dc-voltages, 
if we replace 7z,z+i by its time-average < 7;,;+i >. A junction is called to be in the resistive state if < \i+\{t) >^ 0. 
In the case of one junction in the resistive state there is a finite voltage-drop also in the neighboring junctions due 
to the coupling a. This is shown in Fig. [|. Note that the total voltage is given by 

7= E%i4Etw (19) 
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Figure 3. Real part of a as function of the frequency ui for a d-wave superconductor at two different temperatures 
T. 
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Figure 4. Distribution of the dc-voltage in the neighborhood of one resistive junction as function of the distance 
n and the coupling constant a. 
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Figure 5. Mechanical analog for a stack of coupled Josephson junctions. The angle 7 of a rotator corresponds 
to the phase difference of a josephson junction. One phase is running (= resistive state), the other phases are only 
oscillating (= superconducting state). 



as in the absence of the coupling. 

We want to study in particular the influence of the dynamical coupling between different barriers in the resistive 
state. This will be done both by numerical simulations_and analytical calculations using a Green's function technique 
similar to that used by Takeno in a different context.Ej Details of the calculations will be published elsewhere. 

Let us begin with the discussion of one barrier in the resistive state. It may be helpful to visualize the dynamics 
of such a system by considering the mechanical analog of the RSJ model: The dynamics of each phase difference can 
be described by a pendulum with a constant torque being proportional to the bias current I which is the same in all 
barriers (Fig. |5|). In the absence of the coupling a the pendulum is either rotating, this corresponds to the resistive 
state of the barrier, or has a constant phase, this corresponds to the superconducting state. In the presence of the 
coupling which produces an additional torque, one still can distinguish between these two types of motion: a rotating 
state and a non-rotating vibrating state. In the rotating state there is a running phase, j(t) — 70 + tot + 5~/(t) with a 
finite value of the Josephson frequency u> =< 7 >. In the non-rotating state < f i >= 0, but there are still oscillations. 
Such localized solutions are known in non-linear dynamics as roto-breathers.E£l The running phase in the resistive 
barrier causes finite phase oscillations in the neighboring barriers. The amplitude of these oscillation depends on the 
ratio of the rotation frequency, i.e. the Josephson frequency lu, and the cigcnfrcquencies of the oscillations of the 
coupled system at small amplitudes. The latter are determined by the Josephson plasma frequency uj p and have a 
bandwidth which is proportional to the coupling a. In Fig. ^ we show the result of an analytical calculation of the 
oscillation amplitude for different barriers as a function of the distance from the resistive barrier and as a function of 
the Josephson frequency ui. In this example the plasma frequency is ui p /(2ir) = 0.6 THz and a = 0.2. If the rotation 
frequency is high, u > w p , which is usually the case for intrinsic Josephson systems, the oscillation amplitude in 
the neighboring barriers is small and falls off exponentially with the distance from the resistive barrier. In the case 
of a Josephson frequency u> within the plasma band, the phase-rotation in one barrier leads to long-range plasma 
oscillations in the neighborhood. These are the longitudinal plasma-oscillations considered by several authors which 
in principle can also be excited by longitudinal electric fields in the purely superconducting state. 

From an experimental point of view it is important to know how this effect influences the current-voltage char- 
acteristics of an intrinsic Josephson system, showing a multi-branch structure, where the n-th branch corresponds 
to n resitive barriers. As we have seen above, the total voltage measured for the stack of Josephson junctions is still 
given by the sum of the Josephson frequencies of the resistive barriers. For the total dc- voltage only the values of 
< 7/,;+i > for the resistive barriers contribute. In the case of two or more resistive barriers in the stack it makes a 
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Figure 6. Oscillation amplitude of the phase |7„(a;)| in the different superconducting barriers as function of the 
distance n from the resistive junction and as function of the oscillation frequency uj. 

difference whether these resistive barriers are next neighbors or are separated by one or more non-resistive junctions. 
We have checked this both by numerical simulations and analytical calculations. We find that in the case of two 
resistive barriers next to each other the rotation frequency of both junctions is the same, but (at the same current) is 
slightly higher than in the case of well-separated junctions. For two uncoupled junctions, of course, the total voltage 
is just the double of one resistive junction (the first branch). This is shown in Fig. ^ where we compare the second 
branch of the current voltage characteristics for two neighboring resistive junctions with the total voltage of two well 
separated resistive junctions. For realistic values of a < 1 the effect is very small, therefore we have exaggerated in 
the figure this effect. More pronounced is the difference near the plasma frequency: If we lower the voltage towards 
the upper edge of the plasma band, at least one of the resistive junctions returns to the superconducting state. This 
return point (the minimum of the curves) is different for the two situations: for coupled junctions the return-voltage 
and the dc-current is smaller than for uncoupled resistive junctions. 

Our numerical simulations also show that in the case of two neighboring resistive junctions the two phases rotate 
coherently. This remains true even if the critical currents of the two junctions are slightly different. This phase- 
locking is very important for application of coupled Josephson oscillations for the generation or amplification of 
radiation. 

Our experiments on the intrinsic Josephson systems BSCCO and TBCCO show a very precise additive structure 
of the different branches at least for the lower branches. Deviations from the additivity at higher order branches 
may be attributed to heating effects. This poses an upper limit to the coupling parameter of a < I . An upper limit 
of a is also obtained from the hysteretic return point of the current voltage characteristics from the first resistive 
branch to the superconducting state. A lower limit of this return point is given by the upper edge of the plasma 
band given by uj v \J\ + 2a. From the experimental value for uj p and the return voltage io return one can also estimate 
a < 1. Finally one may calculate a for a 2-dimensional electron gas: for e ss 25 one obtains a = 0.1. 

In this paper we have studied the coupling of Josephson oscillations in different barriers due to charge fluctuations. 
Another mechanism which also leads to such_a coupling is the excitation of c-axis phonons in the stack by Josephson 
oscillations, which we have studied recently£3 The phonons excitated in one resistive barrier induce field oscillations 
also in neighboring barriers and thus support coherence of Josephson oscillations. 



Starting from a microscopic model for the Josephson effect in a stack of superconducting layers coupled by tunneling 
barriers we have derived a set of equations for the gauge-invariant phase difference in different barriers. A simplified 



4. SUMMARY 
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Figure 7. Schematic plot of the current-voltage characteristics for two coupled junctions (solid line) in comparison 
with two uncoupled junctions (dashed line) 

model has then been used to study the non-linear dynamics of the coupled system. In particular, we have studied 
the difference in the current-voltage characteristics of two neighboring resistive junctions in comparison with well 
separated resistive junctions. From the regular additive structure of the current-voltage branches found in experi- 
ments on BSCCO and TBCCO and the return voltage one can derive an upper limit of the coupling constant of the 
order of a < 1. 
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